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[i]  We  investigated  environmental  forcing  mechanisms  of  phytoplankton  absorption  near 
the  Mississippi  River  delta  using  multiyear  satellite  data.  An  algorithm  for  the 
phytoplankton  absorption  coefficient  ( aph )  was  developed  from  in  situ  measurements  and 
applied  to  ocean  color  imagery.  We  employed  a  suite  of  chemical  and  physical  forcing 
variables,  including  surface  currents.  For  satellite-derived  aph  time  series  (2002-2004), 
correlation  and  stepwise  regression  analyses  revealed  the  most  important  forcing  variables 
of  aph  on  the  Louisiana  shelf.  Areally,  Mississippi  River  discharge  and  nitrate 
concentration  ([N03])  were  the  two  most  important  predictors  of  aph  over  the  hypoxic 
zone  (defined  by  its  maximum  extent).  River  discharge  was  important  in  a  band 
stretching  from  the  Mississippi  River  delta  to  the  Louisiana-Texas  border.  Riverine  [N03] 
and  wind  magnitude  best  predicted  aph  in  nearshore  waters,  and  solar  radiation  and 
SST  were  most  important  farther  offshore  over  the  hypoxic  zone,  indicating  upwelled 
nutrient  sources  to  phytoplankton.  A  multiple  linear  regression  model  performed  well  in 
resolving  seasonal  and  interannual  aph  variability  in  model  development  years  (2002- 
2004)  (mean  error  of  1 8%,  over  all  pixels  and  months)  and  in  predicting  shelf-wide  aph 
patterns  in  2005  (mean  error  of  32%).  Our  results  strongly  suggest  that  in  recent 
years,  stratification  and  vertical  mixing,  in  addition  to  riverine  [N03],  play  a  primary  role 
in  regulating  phytoplankton  biomass  over  the  hypoxic  zone.  As  well,  a  springtime  model 
experiment  showed  that  aph  over  the  hypoxic  zone  can  differ  by  an  average  absolute 
37%  from  its  average  scenario  owing  to  changes  solely  in  environmental  variables  other 
than  N03  flux. 

Citation:  Green,  R.  E.,  and  R.  W.  Gould  Jr.  (2008),  A  predictive  model  for  satellite-derived  phytoplankton  absorption  over  the 
Louisiana  shelf  hypoxic  zone:  Effects  of  nutrients  and  physical  forcing,  J.  Geophys.  Res.,  113 ,  C06005,  doi:  10.1 029/2007 JC004594. 


1.  Introduction 

[2]  The  Louisiana-Texas  continental  shelf  is  a  highly 
dynamic  river-dominated  margin  and  the  site  of  the  2nd 
largest  zone  of  coastal  hypoxia  in  the  world.  Increases  in 
nitrate  (N03)  loading  to  the  Mississippi  River  watershed 
during  the  last  50  years  are  considered  responsible  for  an 
increase  in  summertime  hypoxia  in  Louisiana-Texas  shelf 
bottom  waters  [Rabalais  et  al. ,  2002].  There  is  currently  a 
national  mandate  to  decrease  the  size  of  the  hypoxic  zone  to 
5000  km2  by  2015,  mostly  by  a  30%  reduction  in  annual 
nitrogen  flux  into  the  Gulf  of  Mexico  [Mississippi  River/ 
Gulf  of  Mexico  Watershed  Nutrient  Task  Force ,  2001]. 
Hypoxic  bottom  waters  (dissolved  oxygen  concentrations 
<2  mg  L-1)  are  a  concern  primarily  because  of  the  reduced 
catch  of  demersal  fish,  shrimp,  and  crabs  in  these  waters 
[e.g.,  Leming  and  Stuntz ,  1984],  which  otherwise  represent 
one  of  the  United  States’  most  productive  fisheries.  High 
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concentrations  of  dissolved  NO3  from  the  Mississippi  River 
(>100  /iM  l_l  [Dagg  and  Whitledge ,  1991])  support  sig¬ 
nificant  phytoplankton  growth  in  the  river  plume  [Lohrenz 
et  al. ,  1990,  1999].  When  this  primary  production  sinks 
from  surface  to  bottom  waters,  it  is  consumed  by  aerobic 
bacteria  and  is  considered  an  important  factor  leading  to 
hypoxic  conditions.  Given  the  large  areal  extent  of  regional 
hypoxia  (22,000  km2  in  midsummer  2002  [Rabalais  and 
Turner ,  2006]),  a  host  of  environmental  factors  determine 
primary  production  in  surface  waters  over  the  hypoxic  zone 
[e.g.,  Lohrenz  et  al .,  1999].  Assessing  the  mechanisms 
responsible  for  variability  in  phytoplankton  growth  on  the 
shelf  is  requisite  to  a  better  understanding  of  hypoxia 
formation. 

[3]  On  the  basis  of  findings  from  in  situ  measurements, 
the  interaction  of  circulation,  light,  and  nutrient  and  fresh¬ 
water  inputs  all  act  to  regulate  phytoplankton  biomass  and 
production  on  the  Louisiana  shelf.  Seasonal  variability  in 
circulation  and  freshwater  discharge  from  the  Mississippi 
and  Atchafalaya  Rivers  strongly  influence  the  distribution 
of  properties  on  the  shelf.  Maximum  river  discharge  gener¬ 
ally  occurs  in  the  spring  (March-May)  [Dinnel  and 
Wiseman ,  1 986],  and  maximum  retention  of  biomass  and 
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nutrients  on  the  shelf  due  to  favorable  winds  and  currents 
occurs  from  Fall  through  Spring  [Cochrane  and  Kelly , 
1986].  Conditions  of  both  high  discharge  and  favorable 
circulation  have  been  associated  with  higher  phytoplankton 
production  and  biomass  on  the  shelf  [e.g.,  Chen  et  al, 
2000].  As  well,  at  a  midshelf  location  between  the  Mis¬ 
sissippi  and  Atchafalaya  River  mouths,  primary  production 
inferred  from  dissolved  oxygen  concentrations  was  posi¬ 
tively  correlated  with  river  discharge  (lagged  by  2  months) 
over  a  six  year  period  [Justic  et  al. ,  1993].  Linkages 
between  high  concentrations  of  riverbome  nutrients  and 
productivity  are  the  basis  for  management  scenarios  for 
regulating  the  size  of  summertime  hypoxia.  On  the  basis  of 
measurements  from  several  years,  positive  correlations  were 
observed  between  primary  production  and  NO3  +  nitrite 
(N02)  concentrations  near  the  Mississippi  River  delta 
[Lohrenz  et  al. ,  1997].  Other  biological  and  physical  factors 
which  cause  variability  in  phytoplankton  growth  on  the 
shelf  include  mixed  layer  irradiance,  temperature,  and 
salinity  [Lohrenz  et  al .,  1990,  1997,  1999;  Chen  et  al ., 
2000],  as  well  as  phytoplankton  species  composition 
[Redalje  et  al.,  1994]  and  grazing  [Fahnenstiel  et  al ., 
1995].  While  field  measurements  provide  important  insights 
into  oceanographic  processes,  they  are  limited  by  the  spatial 
and  temporal  scales  of  variability  that  they  can  capture. 

[4]  Satellite  remote  sensing  has  proved  a  powerful  tool 
for  studying  oceanographic  processes  over  large  areal 
domains  at  relatively  high  spatial  and  temporal  resolution. 
Imagery  from  the  Coastal  Zone  Color  Scanner  (CZCS) 
ocean  color  satellite  (1978-1985)  provided  the  first  clima¬ 
tological  time  series  of  phytoplankton  concentration  in  the 
Gulf  of  Mexico  [Muller- Karger  et  al .,  1991],  showing  high 
pigment  concentrations  during  the  winter  and  low  concen¬ 
trations  during  the  summer.  Seaward  of  the  shelf,  results 
from  a  combination  of  remote  sensing  observations  and 
one-dimensional  ecosystem  modeling  suggested  that  mixed 
layer  depth  was  the  single  most  important  factor  controlling 
seasonal  pigment  concentrations  [Walsh  et  al,  1989; 
Muller-Karger  et  al .,  1991].  Further  inferences  about  phy¬ 
toplankton  in  the  open  Gulf  of  Mexico  were  made  using 
Sea-viewing  Wide  Field-of-view  Sensor  (SeaWiFS)  ocean 
color  data  and  a  three-dimensional  physical  model  [Jolliff 
et  al .,  2008].  In  addition  to  mixed  layer  depth,  this  analysis 
showed  that  mesoscale  circulation  significantly  modulated 
bio-optical  fields  in  the  Gulf  of  Mexico,  including  phyto¬ 
plankton  pigments.  Surface  distributions  of  suspended  par¬ 
ticulate  matter  (SPM)  have  been  studied  in  more  coastal 
waters  of  the  Gulf  using  satellite  data.  Analysis  of  Advanced 
Very  High  Resolution  Radiometer  (AVHRR)  imagery  (1989— 

1 993)  showed  that  river  discharge  and  wind-forcing  were  the 
main  factors  driving  seasonal  variability  in  SPM  distribu¬ 
tions  near  the  Mississippi  River  delta  [  Walker ,  1996].  More 
recently,  correlation  analyses  of  SeaWiFS  imagery  (1997- 
2000)  identified  two  distinct  regions  of  SPM  variability  near 
the  delta  [Salisbury  et  al .,  2004].  Wind-driven  variability  of 
SPM  was  important  in  shallow  areas,  and  this  region  was 
spatially  separated  from  areas  of  high  discharge-related 
variability. 

[5]  In  recent  years,  significant  effort  has  been  devoted  to 
the  development  of  satellite  ocean  color  algorithms  that  are 
valid  for  retrieving  phytoplankton  biomass  in  coastal 
regions.  In  Case  2  optical  waters,  the  accurate  retrieval  of 


phytoplankton  biomass  from  satellite  imagery  is  complicated 
by  the  presence  of  multiple  seawater  constituents  that  affect 
remote  sensing  reflectance  (Rrs)-  High  concentrations  of 
dissolved  organic  material  and  suspended  sediments  can 
result  in  the  overestimate  of  chlorophyll  concentration  from 
satellites  using  empirical  algorithms  [e.g.,  Darecki  and 
Stramski ,  2004].  However,  the  use  of  in  situ  observations 
to  regionally  correct  bio-optical  algorithms  can  result  in 
significantly  improved  satellite  retrievals.  For  example, 
several  studies  have  validated  satellite  retrievals  with  in 
situ  data  in  coastal  regions  and  have  used  the  resulting 
chlorophyll  algorithms  to  describe  variability  in  phyto¬ 
plankton  biomass,  including  in  Chesapeake  Bay  [Harding 
et  al .,  2005]  and  Massachusetts  Bay  [Hyde  et  al .,  2007].  As 
well,  accurate  satellite  retrievals  of  phytoplankton  biomass 
in  coastal  waters  is  complicated  by  difficulties  with  atmo¬ 
spheric  correction,  which  often  result  in  satellite  retrievals 
that  are  too  low  [e.g.,  Siegel  et  al .,  2000].  One  successful 
technique  for  improving  atmospheric  correction  in  coastal 
regions  applies  an  iterative  method  in  which  NIR  water- 
leaving  radiance  is  estimated  from  water  leaving  radiance  at 
670  nm  [Amone  et  al .,  1998;  Stumpf  et  al .,  2003];  this  NIR 
scheme  improves  retrievals  of  water-leaving  radiances  in 
coastal  areas,  has  been  implemented  in  the  most  recent  (5th) 
global  reprocessing  effort  of  SeaWiFS  at  NASA,  and  was 
used  in  our  image  processing.  As  well,  a  new  absorbing 
aerosol  correction  has  been  developed  (which  we  have  also 
applied  during  our  image  processing)  that  further  improves 
satellite  retrievals  of  phytoplankton  biomass  in  Case  2  waters 
[Ransibrahmanakul  and  Stumpf,  2006].  Taken  together, 
these  advances  in  satellite  image  processing  provide  the 
unprecedented  ability  to  better  understand  phytoplankton 
dynamics  in  the  coastal  zone. 

[6]  In  the  present  study,  our  goal  was  to  determine  the 
main  factors  driving  variability  in  satellite-derived  phyto¬ 
plankton  absorption  near  the  Mississippi  River  delta.  Pre¬ 
vious  satellite-based  time  series  analyses  of  phytoplankton 
biomass  have  focused  mainly  on  waters  seaward  of  the  shelf 
in  the  Gulf  of  Mexico.  We  use  extensive  cruise  data  from 
shelf  waters  in  the  northern  Gulf  of  Mexico  to  extend 
our  understanding  of  phytoplankton  dynamics  to  river- 
dominated  Louisiana  shelf  waters.  For  a  3-year  time  period 
(2002-2004),  we  analyzed  time  series  of  SeaWiFS-derived 
phytoplankton  absorption  coefficients  (aph)  and  a  suite  of 
potential  environmental  forcing  mechanisms.  We  prefer  aph, 
rather  than  chlorophyll  concentration,  as  a  proxy  for  phy¬ 
toplankton  biomass  for  two  main  reasons:  (1)  use  of  aph 
avoids  methodological  controversy  in  how  chlorophyll  is 
measured  (i.e.,  fluorometric  versus  HPLC),  and  (2)  as  an 
optical  property,  aph  should  be  more  directly  related  to  Rrs. 
In  reality,  aph  is  primarily  a  function  of  chlorophyll  a  and 
pheopigments,  but  also  reflects  physiological  states  related 
to  photoadaptation  [Neori  et  al.,  1984]  and  nutrient  stress 
[Heath  et  al.,  1990].  First,  we  used  correlation  and  stepwise 
regression  analyses  to  determine  the  most  important  factors 
driving  aph  variability  and  their  spatial  patterns  of  influence. 
Second,  a  multiple  linear  regression  model  was  developed 
for  aph,  and  the  prediction  capabilities  of  this  model  were 
determined  in  application  to  SeaWiFS  imagery  for  2005. 
With  relevance  to  mandated  decreases  in  NO3  loading  to  the 
Mississippi  River  watershed,  we  also  applied  the  model  to 
calculating  variability  in  shelf  phytoplankton  distributions 
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Figure  1.  Map  of  study  region  in  the  northern  Gulf  of  Mexico,  including  the  Louisiana,  Mississippi, 
and  western  Alabama  coastlines.  Locations  of  sample  sites  for  all  seven  field  experiments  are  indicated 
(crosses).  Boxes  A  and  B  indicate  regions  in  the  Mississippi  and  Atchafalaya  River  plumes,  respectively, 
where  example  average  time  series  were  calculated.  The  box  at  BURL1  indicates  location  of  buoy  wind 
measurements  at  Southwest  Pass,  which  is  the  main  outflow  of  the  Mississippi  River  through  the  delta. 
Also,  a  site  of  long-term  monitoring  for  hypoxia  off  Terrebonne  Bay  is  indicated  (station  C6,  star). 


due  to  environmental  variables  other  than  N03  flux.  This 
analysis  should  help  explain  observed  changes  in  biomass 
that  can’t  be  accounted  for  by  N03-flux  variability  alone. 

2.  Methods 

2.1.  Cruise  Measurements 

[7]  Seven  field  experiments  were  conducted  in  northern 
Gulf  of  Mexico  coastal  waters  of  Louisiana,  Mississippi  and 
Alabama  (Figure  1).  Measurements  were  obtained  during 
two  cruises  on  the  Louisiana  shelf  during  27  April  to  1  May 
2006  aboard  the  R/V  Pelican  and  6-11  September  2006 
aboard  the  OSV  Bold .  Surface  water  samples  were  collected 
using  a  conductivity-temperature-depth  profiler/rosette  sys¬ 
tem  equipped  with  sampling  bottles.  All  samples  were 
processed  onboard  ship  within  2-3  h  of  sample  collection. 
As  well,  two  nearshore  samplings  in  shallower  waters  of  the 
Louisiana  shelf  were  conducted  as  transects  out  of  Vermil¬ 
ion,  Terrebonne,  and  Barataria  Bays  (17-19  May  and  17- 
1 9  July  2007).  Three  field  experiments  were  also  conducted 
off  the  Mississippi/Alabama  coasts  during  20-24  May 
2002,  13  December  2005,  and  6  February  2007.  During 
these  experiments,  samples  were  collected  on  day  trips  and 
processed  that  evening  or  the  next  day  onshore.  Sample 
sites  were  generally  located  inside  the  50  m  isobath,  with 
the  majority  of  stations  within  the  15  m  isobath  (Figure  1). 

[8]  Absorption  coefficients  for  particulate  material  (ap)  in 
surface  water  samples  were  measured  spectrophotometri- 
cally.  Particulate  material  was  collected  on  GF/F  filters 
(nominal  pore  size  of  0.7  /im),  and  ap  was  measured  using 
an  Analytical  Spectral  Devices  (ASD)  fiber  optic-based 
spectrophotometer.  Subsequent  to  the  initial  optical  density 
measurements,  filters  were  extracted  for  20  min  in  hot 
methanol  and  reanalyzed  to  determine  the  residual  particu¬ 


late  absorption  by  detritus  ( ad )  [Kishino  et  al. ,  1985].  The 
absorption  coefficient  due  to  methanol-extractable  phyto¬ 
plankton  pigments  ( aph )  was  estimated  by  difference  be¬ 
tween  the  initial  and  postextraction  measurements  ( ap  —  ad). 
All  measurements  were  relative  to  a  blank  filter  saturated 
with  milli-Q  water,  and  the  optical  density  (OD)  at  850  nm 
for  the  blank  was  subtracted  from  each  sample  spectrum.  A 
path  length  amplification  (0)  factor  was  used  to  correct 
optical  density  measurements  for  particle  concentration  and 
filter  effects.  The  0  factor  was  determined  specifically  for 
our  ASD  spectrophotometer  during  an  experiment  at 
Bigelow  laboratory  [Mitchell  et  al ,  2000],  and  was  param¬ 
eterized  as 

0(X)  =  (0.258+  1.153  x  ODf(X))~\  (I) 

where  ODf  is  the  optical  density  of  the  filter  pad.  For  field 
samples,  duplicate  measurements  were  generally  not 
collected.  However,  for  some  coastal  samples,  heterogene¬ 
ities  were  seen  in  the  particulate  material  on  the  filter  pad, 
and  in  these  cases  a  second  measurement  was  taken 
following  rotation  of  the  pad  on  the  sample  holder.  A  set 
of  122  phytoplankton  absorption  spectra  were  obtained  in 
all. 

2.2.  Remote  Sensing  Imagery 

[9]  Local  Area  Coverage  (LAC)  SeaWiFS  data  of  ~1  km 
resolution  were  collected,  processed,  and  archived  for  the 
Gulf  of  Mexico  region  (80°W-98°W;  17°N-31°N). 
Imagery  was  acquired  using  a  SeaSpace  TeraScan  satellite 
receiving  system  and  processed  with  the  Naval  Research 
Laboratory’s  Automated  Processing  System  (APS) 
[Martinolichy  2005].  APS  Version  3.4  utilized  atmospheric 
correction  algorithms  proscribed  by  NASA’s  fifth  SeaWiFS 
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Table  1.  Overview  of  Chemical  and  Physical  Forcing  Variables  (2002-2004) 


Variable 

Range 

Units 

Source 

Average  |r|“ 

In  Modclb 

%  Pixels^ 

River  Discharge  (0  and  1  month  lag) 

6.2  x  103-2.8  X  104 

m3  s  1 

USACE 

0.26,  0.27 

X,  X 

12,9 

River  [SPM] 

71-298 

mg  L  1 

USGS 

0.20 

X 

5 

River  [N03]  (0  and  1  month  lag) 

0.66-2.4 

mg  L  1 

USGS 

0.28,  0.23 

X,  X 

14.  6 

River  [PO4] 

0.030-0.091 

mg  L  1 

USGS 

0.18 

River  [SiOi] 

4.3 -8.0 

mg  L  1 

USGS 

0.16 

X 

6 

Wind  speed 

3. 7-8.2 

m  s  1 

NDBC 

0.26 

X 

8 

Wind  u  speed 

-4.4- 1.6 

m  s  1 

NDBC 

0.18 

Wind  v  speed 

-3.3-4.2 

m  s  1 

NDBC 

0.21 

Wave  height 

0.47-1.6 

m 

LSU  WAVCIS 

0.15 

X 

4 

Water  height 

0.13-0.52 

m 

LUMCON 

0.22 

X 

7 

Solar  radiation 

230-630 

/xEin  m  2  s  1 

LUMCON 

0.27 

X 

8 

Precipitation 

0.28-18 

mm 

LUMCON 

0.12 

Sea  surface  temperature  (SST) 

12-31 

°C 

IASNFS  model 

0.26 

X 

13 

Current  speed  (SSS) 

0-1.5 

m  s  1 

IASNFS  model 

0.22 

Current  u  speed  (SSU) 

-0.60-1.3 

m  s  1 

IASNFS  model 

0.18 

X 

4 

Current  v  speed  (SSV) 

cn 

1 

© 

T 

m  s  1 

IASNFS  model 

0.16 

X 

4 

“Absolute  correlation  with  aph  time  scries  as  an  average  over  the  entire  hypoxic  zone  (Figure  5). 
bX  denotes  variables  included  in  the  final  multivariate  model. 

cFrom  stepwise  analysis,  the  percent  of  pixels  in  which  each  variable  was  the  most  important  predictor  over  the  hypoxic  zone. 


reprocessing  [Amone  et  al. ,  1998;  Stumpf  et  al. ,  2003].  The 
present  atmospheric  correction  method  used  by  APS 
improves  estimates  of  bio-optical  parameters  in  coastal 
regions  by  applying  an  iterative  technique  in  which  water- 
leaving  radiance  at  765  and  865  nm  is  estimated  from  water 
leaving  radiance  at  670  nm.  An  absorbing  aerosol  correction 
was  also  applied  to  improve  underestimates  of  satellite- 
derived  water  reflectance  [Ransibrahmana/cul  and  Stumpf \ 
2006].  In  our  time  series  analyses,  we  used  monthly 
composite  imagery  for  2002-2005  which  were  temporal 
averages  of  all  valid  pixels  (e.g.,  cloud  free)  in  each  month. 

[10]  A  relationship  was  developed  between  satellite- 
derived  chlorophyll  and  in  situ  phytoplankton  absorption 
measurements.  We  used  a  chlorophyll  product  which  is 
automatically  generated  by  APS,  on  the  basis  of  the 
algorithm  of  Stumpf  et  al.  [2000],  hereafter  referred  to  as 
“Stumpf  chlorophyll.”  This  algorithm  is  parameterized 
using  a  polynomial  regression  of  the  log-transformed  ratio 
of  remote  sensing  reflectance,  Rrs,  at  490  and  555  nm  (rat). 
First,  an  initial  chlorophyll  value  (Q)  in  pg  L~l  is  calcu¬ 
lated  as 

C,  =  io-25xro/.  (2) 

The  following  algorithm  is  then  applied  which  was 
developed  specifically  for  coastal  chlorophyll  concentra¬ 
tions.  This  algorithm  is  a  weighted  transform  which 
blends  to  the  global  OC2v2  chlorophyll  algorithm  in  low- 
chlorophyll  (offshore)  waters: 

if  C,  >  0.5,  then  C  =  C„  (3) 


if  0.1  <  C,  <  0.5,  then 

c=  lOl*0®10  (G)x [|o8io  (G)+1.0]+k>g,0  (CooJxOogio  (0.5) -log, 0  (Ci))}/  l°gio  (5-0)^ 
if  C,  <0.1,  thenC  =  Con  =  10lfl,+rt,,xlfl2+n,/x(as+r,,/x^)l>  -I .a$. 

The  coefficients  a]  to  a5  are  from  the  global  chlorophyll 
algorithm  (OC2v2)  and  equal  0.2974,  -2.2429,  0.8358, 


-0.0077,  and  -0.0929  [Maritorena  and  O'Reilly ,  2000]. 
It  was  not  our  intention  to  validate  the  Stumpf  chlorophyll 
product  for  our  study  area,  but  rather  to  use  satellite-derived 
values  of  this  product  to  determine  a  relationship  with  our 
measurements  of  phytoplankton  absorption  at  443  nm. 
Satellite  matchups  with  sample  locations  were  taken  from 
either  the  same  day  or,  if  not  available,  a  day  either  side  of 
sampling.  Our  derived  relationship  between  satellite-derived 
Stumpf  chlorophyll  and  measured  aph  (443)  was  then  applied 
to  all  satellite  imagery  for  the  time  series  analyses  to  obtain 
satellite  aph  (443)  imagery. 

2.3.  Time  Series  Analyses 

[n]  We  compiled  time  series  of  numerous  chemical  and 
physical  forcing  variables  near  the  Mississippi  River  delta 
(Table  1).  The  following  environmental  variables  were 
chosen  as  those  that  could  plausibly  affect  phytoplankton 
growth  in  surface  waters.  We  used  riverine  discharge  and 
N03  concentration  ([NO3])  as  separate  variables  in  our 
analysis  to  separate  the  physical  effects  of  stratification 
from  biological  effects  of  eutrophication.  Daily  measure¬ 
ments  of  water  discharge  for  the  Mississippi  River  were 
obtained  from  the  United  States  Army  Corps  of  Engineers 
(US ACE)  for  Tarbert  Landing,  Louisiana.  Outflow  from  the 
Mississippi  and  Atchafalaya  Rivers  is  highly  in  phase, 
reflecting  efforts  of  the  USACE  to  manage  the  discharge 
of  the  Mississippi  River  system  so  that  ~2/3  of  the 
discharge  outflows  via  the  bird-foot  delta  and  the  other 
1/3  through  the  Atchafalaya  (Figure  1).  Concentrations  of 
riverine  NO3  ([N03]),  orthophosphate  ([PO4]),  silicate 
([Si03]),  and  SPM  ([SPM])  were  obtained  from  the  United 
States  Geological  Survey  (USGS)  for  St.  Francisville  and 
Baton  Rouge,  Louisiana  (Figure  1);  data  were  generally 
available  bimonthly  for  each  station.  Daily  values  were 
interpolated  from  the  USGS  data  and  then  monthly  means 
were  calculated.  Two  additional  variables  of  riverine  [N03] 
and  discharge  with  a  one  month  lag  were  also  included  in 
our  analysis.  Hourly  measurements  of  wind  speed  and 
direction  were  obtained  from  the  National  Data  Buoy 
Center  (NDBC)  for  a  station  near  the  main  Mississippi 
River  outflow  at  Southwest  Pass  (BURL1;  Figure  1).  We 
analyzed  wave  height  data  from  Louisiana  State  Univer- 
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Figure  2.  Algorithm  for  phytoplankton  absorption  (aph)  at 
443  nm  as  a  function  of  satellite-derived  chlorophyll. 
Stumpf  chlorophyll  was  calculated  from  satellite-derived 
/?rA.(490)//?r5(555)  as  described  in  section  2.  Phytoplankton 
absorption  values  were  determined  from  filter  pad  absorp¬ 
tion  measurements  collected  during  seven  cruises  in  the 
northern  Gulf  of  Mexico.  Note  that  the  coefficient  of 
determination  ( r 2  value)  explains  variance  in  log- 
transformed  aph. 


sity’s  (LSU’s)  Wave-Current-Surge  Information  System 
(WAVCIS)  for  Station  5  (-90°32',  29°3.2'),  which  is 
located  south  of  Terrebonne  Bay.  Several  variables  were 
obtained  from  the  Louisiana  Universities  Marine  Consor¬ 
tium  (LUMCON)  weather  station,  including  solar  radiation, 
precipitation,  and  water  height. 

[12]  A  physical  circulation  model  provided  sea  surface 
temperature  (SST),  current  magnitude  (SSS),  and  current 
direction  (SSU,  SSV).  We  used  the  IntraAmericas  Sea  Now¬ 
cast/Forecast  System  (IASNFS)  [Ko  et  al. ,  2003]  which 
consists  of  a  1/24°  (5-6  km),  41 -level,  data-assimilating 
ocean  model  based  on  the  Navy  Coastal  Ocean  Model 
(NCOM)  [Martin,  2000].  The  IASNFS  includes  the  Gulf  of 
Mexico,  the  Caribbean  Sea,  the  Straits  of  Florida,  and  part  of 
the  western  North  Atlantic  Ocean.  The  ocean  model  contin¬ 
uously  assimilates  three-dimensional  ocean  temperature/ 
salinity  analyses  produced  by  a  statistical  data  analysis 
model,  the  Modular  Ocean  Data  Assimilation  System 
(MODAS)  [ Carnes  etal. ,  1996].  Real-time  satellite  altimeter 
sea  surface  height  anomaly  data,  as  well  as  AVHRR  SST  data 
are  assimilated  into  the  model.  Surface  forcing  is  provided  by 
surface  wind  stress,  heat  fluxes,  solar  radiation,  air  pressure, 
and  river  discharge  (monthly  climatology)  from  53  rivers. 
Rivers  are  prescribed  by  river  temperature  (monthly  clima¬ 
tology),  salinity  (fresh  water),  and  discharge  as  a  function  of 
depth.  The  open  boundary  conditions  are  provided  by  baro- 
tropic  tidal  forcing  and  the  Navy’s  global  NCOM  model  [Ko 
et  al .,  2007].  We  define  currents  and  winds  in  the  u  (zonal) 
and  v  (meridional)  directions,  where  positive  u  is  to  the  east 
and  positive  v  is  to  the  north.  For  our  analysis,  we  interpo¬ 
lated  SST  and  currents  to  the  higher  resolution  of  the 
SeaWiFS  imagery  in  a  region  encompassing  the  Louisiana 
shelf  (latitudes  28.4°N  to  30°N,  longitudes  89°W  to  94°W). 
Mean  monthly  values  were  calculated  for  all  environmental 
variables. 


[13]  Correlation  and  regression  analyses  were  performed 
to  analyze  the  impacts  of  the  different  forcing  variables  on 
driving  mean  monthly  optical  variability.  In  the  case  of 
point  measurements,  such  as  river  discharge  and  wind 
speed,  the  same  physical  time  series  was  assumed  at  all 
image  pixels.  In  contrast,  a  unique  physical  time  series  at 
each  image  pixel  was  used  for  IASNFS  model  results. 
Spikes  were  removed  from  satellite-derived  aph( 443)  time 
series  by  discarding  points  that  were  more  than  five  times 
the  mean  at  each  pixel  over  the  3  year  study  period  (2002- 
2004).  As  a  first  step,  correlation  analyses  were  performed 
between  mean  monthly  satellite-derived  aph( 443)  time 
series  for  2002  2004  and  each  forcing  variable.  Next,  we 
used  stepwise  multiple  linear  regression  analysis  to  deter¬ 
mine  the  importance  of  the  different  variables  in  determin¬ 
ing  aph  variability.  A  forward  stepwise  procedure  was  used 
in  which  the  most  statistically  significant  term  (the  one  with 
the  lowest  p  value  for  p  <  0.05)  is  added  into  the  model  at 
each  step,  until  there  are  no  statistically  significant  terms  left 
to  include.  A  common  problem  with  multiple  regression 
analysis  is  multicollinearity  of  the  input  variables,  a  point 
which  we  will  further  discuss.  On  the  basis  of  the  stepwise 
analysis,  several  variables  were  identified  as  least  important 
and  were  removed  from  our  final  multivariate  model. 

[14]  Using  the  2002-2004  data  sets,  a  multiple  linear 
regression  model  for  aph  was  developed  at  each  pixel,  such 
that  forcing  variables  had  different  model  weightings 
depending  on  location.  The  resulting  model  was  then 
applied  to  predicting  aph  for  several  different  scenarios. 
Mean  monthly  environmental  variables  for  2005  were  input 
to  the  model  to  determine  how  well  temporal  and  spatial 
variability  in  aph  could  be  predicted  for  a  year  that  was  not 
included  in  model  development.  As  well,  we  performed  an 
experiment  with  the  model  to  see  how  much  aph  could 
change  in  spring  (April— June)  2002-2004  on  the  basis  of 
natural  variability  in  all  environmental  variables,  except 
nitrate  flux  (riverine  [NO3]  x  discharge).  For  this  experi¬ 
ment,  a  maximum  and  minimum  aph  scenario  were  modeled 
using  the  measured  limits  of  all  springtime  environmental 
variables,  while  both  riverine  [NO3]  and  discharge  were 
held  constant  at  their  average  values  for  this  time  period. 
Decreases  in  nitrogen  flux  from  the  Mississippi  River  to  the 
Gulf  of  Mexico  have  been  mandated  in  order  to  reduce  the 
size  of  summertime  hypoxia  [Mississippi  River/Gulf  of 
Mexico  Watershed  Nutrient  Task  Force ,  2001].  As  riverine 
N03  flux  is  reduced,  certainly  other  environmental  factors 
will  cause  some  variability  in  phytoplankton  biomass  on  the 
shelf  and  the  size  of  hypoxia.  Thus,  we  wanted  to  explore 
how  much  these  other  sources  of  variability  could  affect  ap^ 
The  experiment  was  performed  for  the  spring,  a  period 
during  which  in  situ  phytoplankton  growth  in  surface  waters 
is  considered  important  in  the  development  of  regional 
hypoxia.  We  performed  all  time  series  data  analyses  using 
the  MATLAB  software  package  (The  Math  Works). 

3.  Results 

3.1.  Algorithm  Development 

[15]  We  developed  a  regional  algorithm  for  phytoplank¬ 
ton  absorption  at  443  nm  based  on  satellite-derived  remote 
sensing  reflectance  for  the  northern  Gulf  of  Mexico 
(Figure  2).  Measurements  of  ^(443)  on  the  Louisiana 
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Figure  3.  For  2002-2005,  time  series  of  (a)  satellite-derived  aph  (443),  (b)  Mississippi  River  discharge, 
(c)  wind  speed,  and  (d)  riverine  [N03].  The  aph  time  series  presented  here  are  an  average  over  Box  A  in 
the  Louisiana  Bight  (Figure  1).  The  first  3  years  (2002-2004)  were  used  to  develop  the  optical  model, 
which  was  subsequently  applied  to  2005  data  to  predict  aph  for  that  year. 


shelf  from  all  four  field  experiments  ranged  from  0.02  to 
2.9  m^1.  Peak  values  were  measured  in  Barataria  Bay 
during  nearshore  sampling  in  May  and  July,  2007.  Of  the 
122  measured  aph( 443)  values,  there  were  only  63  matchups 
with  satellite  imagery  owing  to  clouds  over  the  sample  site. 
There  was  substantial  variability  in  measured  aph( 443)  for 
the  same  value  of  satellite-derived  chlorophyll  concentra¬ 
tion  (Figure  2),  and  this  variability  was  even  within  a  single 
field  experiment.  This  variability  in  ^(443)  averaged 
~2. 5-fold,  similar  to  previous  observations  in  these  waters 
[D'Sa  et  al. ,  2006].  A  good  fit  was  obtained  between 
satellite-derived  Stumpf  chlorophyll  and  log-transformed 
values  of  measured  ^(443)  using  a  power  function  (r2  = 
0.85;  RMSE  =  0.18).  This  relationship  was  used  to  calculate 
the  satellite-derived  ^(443)  product  used  in  the  following 
time  series  analyses. 

3.2.  Physical  Time  Series 

[i6]  For  2002-2005,  phytoplankton  absorption  peaked 
from  late  winter  to  summer  (March— August)  on  the  Louisiana 
shelf,  with  significant  spatial  and  interannual  variability 
observed  in  peak  magnitude  and  timing.  Highest  aph 
values  were  usually  associated  with  the  Mississippi  and 
Atchafalaya  River  outflows,  though  locally  high  values  were 
also  observed  outside  of  Barataria  Bay  and  near  the  Sabine 
River  outflow.  To  exemplify  temporal  trends,  we  analyzed 
aph  in  two  different  regions  of  the  shelf:  the  Louisiana  (LA) 
Bight  in  the  Mississippi  River  plume  (Box  A;  Figure  1)  and 
the  Atchafalaya  River  plume  (Box  B).  Values  of  aph  in 
Boxes  A  and  B  ranged  from  0.08  to  0.30  m-1  and  0.17  to 
0.79  m_1,  such  that  aph  in  the  Atchafalaya  plume  box 
was  ~2  times  higher  than  in  the  Mississippi  plume  box. 
Differences  in  peak  timing  were  also  observed  between  the 
two  regions.  Whereas  aph  in  the  Atchafalaya  plume  always 
peaked  in  July  (not  shown),  maxima  in  the  Mississippi 


plume  were  observed  from  March  through  August,  depend¬ 
ing  on  the  year  (Figure  3a).  As  well,  interannual  and  spatial 
variability  was  observed  in  whether  aph  had  one  or  two 
annual  peaks.  A  notable  exception  to  peak  timing  was 
occasionally  observed  in  the  southernmost  region  of  the 
hypoxic  zone,  south  of  29°N  latitude  (extent  of  hypoxic  zone 
shown  in  Figure  5).  In  this  region,  the  influence  of  offshore 
waters,  most  likely  in  the  form  of  upwelled  nutrients,  occa¬ 
sionally  caused  aph  to  peak  in  winter  (December- January). 

[n]  Most  environmental  forcing  variables  demonstrated 
some  degree  of  seasonal  trend  during  the  period  of  our 
study.  As  expected,  highest  values  were  observed  for  river 
discharge  in  winter  to  spring  (Figure  3b),  for  wind  speed  in 
winter  (Figure  3c),  and  for  riverine  [N03]  in  spring  to 
summer  (Figure  3d).  Springtime  maxima  were  also  gener¬ 
ally  observed  in  solar  radiation,  precipitation,  and  wave 
height,  and  summertime  maxima  were  observed  in  [P04] 
and  water  height  and  for  SST  in  the  Louisiana  Bight. 
Wintertime  peaks  were  observed  in  both  [SPM]  and 
[Si03].  Both  wind  and  current  u  speed  (SSU  in  LA  Bight) 
were  largest  to  the  west  in  winter,  and  wind  and  current  v 
speed  (SSV  in  LA  Bight)  to  the  north  in  spring.  Throughout 
much  of  the  year,  winds  and  Ekman  transport  caused 
circulation  to  flow  westward  over  the  Louisiana  shelf 
(e.g.,  Figure  4)  [Cochrane  and  Kelly ,  1986].  However, 
between  November  and  February,  this  general  pattern  was 
periodically  altered  by  N  or  NW  winds  associated  with  the 
passage  of  continental  air  masses  causing  offshore  circula¬ 
tion  over  the  Louisiana  shelf.  Also,  in  July  offshore  flow 
was  always  observed  owing  to  periods  of  W  winds. 

[is]  A  correlation  matrix  was  computed  for  all  time  series 
(Table  2).  In  the  Mississippi  River  plume  (Figure  1,  Box  A), 
aph  was  significantly  correlated  with  five  forcing  variables, 
including  (in  order  of  importance)  discharge  (r  =  0.59), 
[N03]  (r  =  0.51),  lagged  discharge  (r  =  0.49),  wind  speed 
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Figure  4.  Example  IntraAmericas  Sea  Nowcast/Forecast 
System  (IASNFS)  circulation  model  fields  for  sea  surface 
temperature  (SST),  current  speed  (SSS),  and  current 
direction  as  a  mean  for  March  2002.  This  example  is 
representative  of  mean  westward  flow  over  the  Louisiana 
shelf  during  much  of  the  year. 


(r  =  -0.33),  and  solar  radiation  (r  =  0.33).  Of  the  171 
possible  correlations  between  the  18  forcing  variables,  33% 
were  considered  significant,  with  an  absolute  r  value  >  0.33 
(two-tailed  test,  a  =  0.05).  Particularly  high  correlations 
(|r|  >  0.70)  were  observed  between  [N03]  and  lagged 
discharge,  [P04]  and  SST,  wind  speed  and  SST,  and  SSS 
and  SSU.  In  contrast  to  the  Mississippi  River  plume,  aph  in 
the  Atchafalaya  River  plume  (Figure  1,  Box  B)  was  most 
highly  correlated  with  [N03]  (r  =  0.59),  followed  by  wind 
speed  (r  =  —0.55)  (data  not  shown). 

3.3.  Correlation  Analyses 

[19]  We  first  performed  simple  correlation  analyses  be¬ 
tween  time  series  of  aph  and  each  environmental  forcing 
variable  for  2002-2004.  An  average  correlation  coefficient 
for  each  variable  was  calculated  over  the  hypoxic  zone 


(Table  1  and  Figure  5).  This  zone  is  here  defined  as  the 
region  from  the  coast  to  the  maximum  seaward  extent  of 
hypoxic  waters  in  summer  based  on  measurements  from 
1985  to  2002  [Turner  et  al .,  2006].  Average  correlations 
showed  that  six  variables  were  most  important  in  determin¬ 
ing  aph  variability  (|r|  >  0.26),  including  river  discharge 
with  0  and  1  month  lags,  riverine  [N03]  (0  lag),  wind  speed, 
SST,  and  solar  radiation  (Table  1  and  Figure  5).  Average 
correlations  were  generally  low,  because  values  were  com¬ 
puted  over  the  entire  hypoxic  zone  (Table  1).  However, 
much  higher  values  (up  to  0.85)  were  observed  in  some 
regions  (Figure  5).  In  nearshore  waters,  aph  was  positively 
correlated  with  all  six  most  important  variables,  except  wind 
speed  with  which  it  was  negatively  correlated.  In  offshore 
waters,  but  still  over  the  hypoxic  zone,  aph  was  negatively 
correlated  with  SST  and  solar  radiation.  Thus,  two  different 
modes  of  forcing  were  observed  for  SST  and  solar  radiation 
depending  on  cross-shelf  location. 

[20]  Although  of  less  overall  importance,  the  other  envi¬ 
ronmental  variables  in  our  analysis  were  significantly  cor¬ 
related  (|r|  >  0.33)  with  aph  somewhere  in  our  region  of 
interest.  For  example,  riverine  [N03]  lagged  by  1  month 
roughly  showed  correlation  patterns  (not  shown)  similar  to 
[N03]  (0  lag),  though  of  slightly  lower  magnitude  (Table  1 ). 
Higher  wind  and  current  speeds  toward  the  east  were 
positively  and  significantly  correlated  with  aph  in  the 
Atchafalaya  plume.  This  corresponds  to  our  previous  ob¬ 
servation  of  maximum  aph  occurring  in  the  Atchafalaya 
plume  in  July,  a  month  when  currents  were  always  east¬ 
ward.  Higher  wind  and  current  speeds  toward  the  north 
were  positively  and  significantly  correlated  with  aph  in  the 
Louisiana  Bight,  reflecting  the  increased  impact  of  the 
Mississippi  River  plume  under  these  conditions.  Water 
height  showed  significant,  but  negative,  correlations  with 
aph  in  a  region  south  of  Atchafalaya  and  Terrebonne  Bays 
between  the  coast  and  50  m  isobath.  These  findings  em¬ 
phasize  the  high  degree  of  spatial  variability  observed  in  the 
environmental  factors  which  drive  aph  in  a  region  as  large  as 
the  Louisiana  shelf. 


Table  2.  Correlation  Matrix  for  aph  and  Forcing  Variables  in  the  Louisiana  Bight  (2002-2004)8 


aph 

Dis 

DisL 

SPM 

[NO?] 

[NOi3]L 

[P04] 

[SiOj] 

Wnds 

Wndu 

Wndv 

WavH 

WatH 

Rad 

Prec 

SST 

SSS 

SSU 

SSV 

1.00 

Dis 

0.59 

1.00 

DisL 

0.49 

0.53 

1.00 

SPM 

0.32 

0.51 

0.21 

1.00 

[NO,] 

0.51 

0.48 

0.70 

0.15 

1.00 

[NOi]L 

0.27 

0.18 

0.48 

-0.17 

0.52 

1.00 

[PO4] 

-0.14 

-0.42 

-0.11 

-0.52 

0.14 

0.29 

1.00 

[SiO,] 

0.05 

0.01 

-0.09 

0.01 

0.19 

-0.04 

0.40 

1.00 

Wnds 

-0.33 

-0.09 

-0.42 

0.37 

-0.50 

-0.66 

-0.39 

0.18 

1.00 

Wndu 

0.17 

0.16 

0.17 

-0.15 

0.44 

0.27 

0.06 

0.17 

-0.41 

1.00 

Wndv 

0.19 

0.19 

0.51 

0.00 

0.46 

0.25 

0.18 

-0.20 

-0.58 

0.04 

1.00 

WavH 

0.07 

0.29 

0.08 

0.44 

0.05 

-036 

-0.26 

-0.09 

0.43 

-0.24 

0.19 

1.00 

WatH 

-0.10 

-0.25 

-0.20 

-0.29 

-0.13 

0.20 

0.49 

-0.14 

-033 

-0.13 

0.12 

-0.05 

1.00 

Rad 

033 

0.28 

0.61 

-0.08 

0.51 

0.58 

0.21 

-034 

-0.68 

-0.01 

0.68 

0.02 

0.36 

1.00 

Prec 

-0.20 

-0.06 

-0.09 

-0.25 

0.01 

-0.07 

0.30 

0.08 

-0.20 

0.10 

035 

-0.01 

0.25 

-0.06 

1.00 

SST 

-0.03 

-0.27 

0.10 

-0.57 

0.23 

0.48 

0.83 

0.05 

-0.73 

0.13 

0.50 

-0.25 

0.64 

0.58 

0.34 

1.00 

SSS 

0.19 

0.33 

0.10 

0.29 

-0.07 

0.14 

-0.15 

-0.18 

0.09 

-0.58 

0.12 

0.32 

0.13 

039 

-0.21 

-0.05 

1.00 

SSU 

0.03 

-0.14 

0.16 

-0.23 

0.37 

0.05 

0.20 

0.17 

-0.29 

0.65 

0.15 

-0.26 

-0.08 

-0.11 

0.20 

0.17 

-0.87 

1.00 

SSV 

0.29 

0.35 

0.20 

0.14 

0.26 

0.30 

-0.05 

-0.31 

-0.33 

-0.23 

034 

0.20 

0.12 

0.49 

-0.05 

0.14 

0.55 

-0.28 

1.00 

“Correlations  in  this  tabic  for  apfn  SST,  SSS,  SSU,  and  SSV  arc  averaged  over  Box  A  (Figure  1).  Absolute  values  greater  than  the  critical  value  of  0.33 
(a  =  0.05)  arc  in  bold.  N  =  36.  Variables:  Dis,  Mississippi  River  discharge;  DisL,  1  month  lagged  river  discharge;  [N03]L,  1  month  lagged  [NO3];  Wnds, 
wind  speed;  Wndu,  wind  u  speed;  Wndv,  wind  v  speed;  WavH,  wave  height;  WatH,  water  height;  Rad,  solar  radiation;  Prec,  precipitation. 
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Figure  5.  For  mean  monthly  time  series  (2002  -2004),  absolute  correlation  of  satellite-derived  aph  (443) 
with  (a)  Mississippi  River  discharge,  (b)  river  discharge  with  1  month  lag,  (c)  riverine  [N03],  (d)  wind 
speed,  (e)  solar  radiation,  and  (0  SST.  The  maximum  seaward  extent  of  the  hypoxic  zone  is  indicated 
(dashed  line)  [Turner  et  al.y  2006], 


3.4.  Phytoplankton  Absorption  Model 

[21]  Stepwise  multiple  linear  regression  analysis  was  used 
to  determine  the  importance  of  the  various  environmental 
variables  in  determining  aph  on  the  shelf.  Only  significant 
variables  (p  <  0.05)  were  retained  in  the  model  (see 
section  2).  We  calculated  the  importance  of  each  variable 
based  on  the  number  of  pixels  in  which  it  was  included  in 
the  stepwise  model  over  the  hypoxic  zone.  The  six  most 
important  variables  were  [N03]  (0  lag),  discharge  with  0 
and  1  month  lags,  SST,  solar  radiation,  and  wind  speed:  the 
same  variables  identified  in  the  correlation  analyses.  At 
each  pixel,  we  determined  which  of  these  variables  was  the 
best  predictor  of  aph ,  and  we  observed  coherent  spatial 
patterns  across  the  shelf  (Figure  6).  Discharge  and  [N03] 
(including  lagged  variables)  were  the  most  important 
predictors  of  aph  over  the  hypoxic  zone  and  were  approx¬ 
imately  equal  in  importance.  Discharge  was  spatially  the 
best  predictor  along  a  band  stretching  from  the  Mississippi 
River  delta  to  the  Sabine  River.  Nitrate  concentration  and 
wind  magnitude  tended  to  be  the  best  predictors  of  aph  in 
shallow,  nearshore  waters,  whereas  solar  radiation  and  SST 
were  the  best  predictors  farther  offshore.  Although  [N03] 
appears  to  be  important  in  some  waters  farther  offshore 
(Figure  6),  in  these  regions  the  relationship  between  [N03] 
and  aph  was  generally  not  significant  (p  >  0.05).  On  the 
basis  of  the  stepwise  analysis,  we  also  found  that  certain 
parameters  were  not  important  to  include  in  the  aph  model, 
including  wind  u  and  v  magnitude,  [PO4],  precipitation, 
and  current  speed.  These  five  parameters  were  least  often 
included  in  the  stepwise  model  and  generally  did  not  show 


coherent  spatial  patterns  in  their  model  inclusion;  they 
were  removed  from  the  following  analyses. 

[22]  A  multiple  linear  regression  model  was  developed 
for  phytoplankton  absorption  at  each  pixel  using  the  remain¬ 
ing  environmental  variables  (Table  1).  Thus,  different  model 
weightings  were  assigned  to  forcing  variables  depending  on 
location.  We  then  applied  the  model  to  the  environmental 
data  to  create  monthly  aph  distribution  maps  for  comparison 
with  the  satellite  monthly  composite  aph  imagery.  For  mean 
monthly  time  series  (2002-2004),  the  average  correlation 
coefficient  (r)  between  modeled  and  satellite-derived  aph 
was  0.82  over  our  region  of  interest  with  all  variables 
included.  The  five  least  important  parameters  (wind  u  and 
v  magnitude,  [PO4],  precipitation,  and  current  speed)  were 
removed  to  develop  the  most  parsimonious  model.  The 
resulting  multiple  regression  model  performed  well,  with 
modeled  and  satellite-derived  aph  time  series  that  were 
significantly  correlated  at  all  pixels  and  with  an  average  r 
of  0.75.  The  main  seasonal  and  interannual  features  in  aph 
were  generally  resolved  (e.g.,  Figures  7  and  8).  However, 
there  were  occasional  differences  between  satellite  and 
modeled  aph ,  such  as  in  the  westward  orientation  of  the 
Atchafalaya  River  plume  in  October  2002  (Figures  7g 
and  7h)  and  in  the  seaward  extent  of  the  Mississippi  River 
plume  in  April  2003  (Figures  8c  and  8d).  The  overall 
difference  between  satellite  and  modeled  aph  was  18%, 
averaged  over  all  months  and  pixels. 

[23]  The  model  performed  well  in  resolving  seasonal  and 
interannual  monthly  differences  in  aph  for  the  years  used  in 
model  development.  In  the  Louisiana  Bight  (Box  A, 
Figure  1),  seasonal  maxima  in  satellite-derived  aph  for 
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Figure  6.  For  mean  monthly  time  series  (2002-2004), 
regions  where  each  variable  listed  was  the  most  important  in 
driving  aph( 443).  This  was  determined  from  stepwise 
multiple  linear  regression  analysis  for  the  most  important 
variables  overall:  discharge  (Disch),  wind,  [NO3],  solar 
radiation  (Rad),  and  SST.  In  this  figure,  discharge  and 
[NO3]  include  both  the  0  and  1  month  lagged  variables.  A 
smoothed  image  for  eased  interpretation  was  produced  by 
repeatedly  passing  a  3  x  3  window  over  the  image  and 
assigning  the  most  likely  value  to  the  center  pixel.  The 
seaward  extent  of  the  hypoxic  zone  is  indicated  (dashed 
line),  as  is  a  site  of  long-term  hypoxia  monitoring  (station 
C6,  star). 

2002-2004  were  observed  from  March  to  July  and  minima 
from  August  to  December  (Figure  3a).  These  seasonal 
trends  were  also  seen  in  the  modeled  aph  distributions 
(e.g..  Figure  7).  For  example,  in  2002,  aph  in  the  Mississippi 
and  Atchafalaya  plumes  (Boxes  A  and  B)  were  both  higher 
in  July  than  in  October  (Figures  7e— 7h).  We  adjusted 
values  of  the  different  forcing  variables  to  determine  which 
environmental  variables  were  primarily  responsible  for 
observed  differences  in  aph.  In  the  Mississippi  plume,  lower 
wind  magnitude  in  July  2002  (Figure  3c)  was  one  of  the 
primary  causes  of  higher  aph  than  in  October  2002.  In 
contrast,  in  the  Atchafalaya  plume,  higher  [N03]  in  July 
was  the  main  cause  of  higher  aph  than  in  October.  In 
addition  to  seasonal  patterns,  interannual  monthly  differ¬ 
ences  in  aph  were  also  represented  by  model  results  (e.g.. 
Figure  8).  For  example,  compared  to  2003  and  2004,  spatial 
patterns  in  April  2002  were  characterized  by  aph  that  was 
higher  in  the  Mississippi  plume  (Box  A)  and  lower  in  the 
Atchafalaya  plume  (Box  B)  in  both  satellite  and  modeled 
aph.  In  the  Mississippi  plume,  lower  river  discharge  and 
higher  wind  magnitude  in  both  2003  and  2004  (Figures  3b 
and  3c)  noticeably  decreased  aph  compared  to  2002.  In 
contrast,  in  the  Atchafalaya  plume,  higher  SSV  in  2003  and 
higher  [N03]  in  2004  (Figure  3d)  were  the  primary  reasons 
for  higher  aph  than  in  2002. 

3.5.  Model  Predictions 

[24]  Our  ultimate  goal  was  to  develop  a  model  for 
predicting  phytoplankton  absorption  at  monthly  timescales, 
in  a  year  that  was  independent  of  model  development.  To 
this  end,  we  applied  the  regression  model  to  2005  environ¬ 
mental  data.  Application  of  the  model  to  2005  was  a  good 
test  of  prediction  capabilities,  as  river  discharge  (and  certain 


other  environmental  variables)  was  often  outside  the  ranges 
that  went  into  model  development  (Figures  3b -3d).  These 
instances  often  lead  to  modeled  values  of  aph  which  were 
negative.  To  correct  for  this,  if  an  environmental  variable  in 
2005  was  higher  (or  lower)  than  the  range  for  that  variable 
in  model  development,  we  substituted  it  with  the  maximum 
(or  minimum)  value  for  that  particular  month  from  2002  to 
2004.  The  resulting  model  for  2005  lead  to  a  mean  error 
between  satellite  and  modeled  aph  of  32%  over  all  pixels 
and  months,  which,  as  expected,  was  higher  than  the  1 8% 
error  we  observed  for  2002-2004.  For  comparison,  the 
average  error  between  satellite  imagery  and  aph  climatology 
(1998-2007)  was  also  calculated.  This  error  with  climatol¬ 
ogy  was  slightly  lower  (27%)  than  from  the  model  (32%), 
but  the  model  provides  the  additional  benefit  of  elucidating 
the  factors  driving  optical  variability,  which  cannot  be 
gained  from  the  climatology  alone.  Expected  patterns  in 
shelf-wide  aph  were  generally  reproduced  by  our  model  in 
2005  (e.g.,  Figure  9).  For  example,  the  model  helped  explain 
higher  values  of  satellite-derived  aph  on  the  Louisiana  shelf 
early  in  2005,  compared  to  late  in  the  year.  In  the 
Mississippi  plume  (Box  A),  we  identified  higher  river 
discharge  and  [Si03]  in  January  2005  (e.g..  Figure  3b)  as 
important  drivers  of  the  higher  aph  observed  in  January 
compared  to  December  (Figure  3a).  In  the  Atchafalaya 
River  plume  (Box  B),  higher  lagged  river  discharge  in 
January  2005  was  primarily  responsible  for  higher  aph 
compared  to  December  2005.  As  well,  mean  monthly  time 
series  of  satellite-derived  and  model-predicted  aph  compared 
well  in  2005.  For  example,  both  satellite  and  modeled  aph 
peaked  in  spring  to  early  summer  in  the  Mississippi  River 
plume  (Box  A;  Figure  10a),  in  comparison  to  midsummer 
peaks  in  the  Atchafalaya  River  plume  (Box  B;  Figure  10b). 

[25]  We  also  did  an  experiment  with  the  model  to  deter¬ 
mine  the  impact  environmental  factors  other  than  NO3  flux 
have  on  the  distribution  and  magnitude  of  aph  on  the 
Louisiana  shelf.  Environmental  data  were  changed  within 
the  natural  variability  observed  for  spring  (April— June) 
2002-2004,  while  N03  flux  (riverine  [N03]  x  discharge) 
was  held  constant  at  its  average  value  for  that  period. 
Previously,  we  considered  the  separate  effects  of  riverine 
discharge  and  nutrient  concentration.  However,  in  this 
experiment,  we  consider  the  two  together  as  NO3  flux, 
because  this  is  the  quantity  specified  by  mandated  nutrient 
reductions  [Mississippi  River/Gulf  of  Mexico  Watershed 
Nutrient  Task  Force ,  2001].  A  springtime  scenario  using 
average  environmental  data  (Figure  11a)  was  compared  to 
maximum  and  minimum  aph  scenarios  (Figures  lib 
and  11c).  These  scenarios  were  modeled  using  the  limits 
of  all  springtime  environmental  variables,  except  NO3  flux. 
On  the  basis  of  our  model  for  2002-2004,  the  results 
suggested  that  springtime  aph  over  the  hypoxic  zone  could 
vary  by  an  average  absolute  difference  of  37%  from  its  ave¬ 
rage  aph  scenario  owing  to  changes  in  environmental  varia¬ 
bles  other  than  NO3  flux.  The  three  most  important  non  flux 
variables  causing  variability  in  the  spring  were  [Si03],  SST, 
and  solar  radiation,  which  contributed  25%  of  the  37% 
change  in  aph.  Much  higher  variations  in  aph  were  observed 
in  specific  regions  of  the  shelf.  For  example,  in  some  regions 
of  the  Louisiana  Bight  (Figure  1),  aph  varied  by  an  average 
absolute  difference  of  ~100%  from  its  average  aph  scenario. 
In  a  second  experiment,  riverine  [NO3]  (versus  flux)  was  held 
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Figure  7.  For  2002,  multiple  linear  regression  model  performance  comparing  satellite-derived  aph  (left- 
hand  column)  and  modeled  aph  (right-hand  column)  in  an  example  from  each  season.  On  each  plot,  the 
location  of  boxes  A  and  B  (Figure  1)  in  the  Mississippi  and  Atchafalaya  River  plumes,  respectively,  is 
indicated.  White  pixels  in  SeaWiFS  imagery  represent  algorithm  failures  (i.e.,  chlorophyll  retrieval 
failure,  so  subsequent  aph  calculation  also  failed). 


constant,  and  average  springtime  aph  over  the  hypoxic  zone 
varied  by  an  average  absolute  difference  of  42%  from  its 
average  aph  scenario  owing  to  changes  in  all  other  environ¬ 
mental  variables  (including  discharge  in  this  case). 

4.  Discussion 

[26]  Our  understanding  of  phytoplankton  dynamics  from 
the  regression  model  is  only  as  good  as  the  accuracy  of 
satellite  retrievals  in  coastal  regions.  In  coastal  waters,  the 
accurate  retrieval  of  phytoplankton  absorption  is  complicat¬ 
ed  by  high  concentrations  of  dissolved  organic  material  and 
suspended  sediments,  components  which  may  sometimes 
co vary  with  phytoplankton  concentrations.  Satellite  algo¬ 
rithm  errors  are  also  caused  by  differences  in  measurement 
times  between  in  situ  collection  and  satellite  overpasses 
[e.g.,  Chang  and  Gould ,  2006].  For  example,  our  absorption 
measurements  were  made  at  various  times  of  the  day, 
whereas  the  SeaWiFS  overpass  occurs  near  local  noon.  As 
well,  if  the  measurement  day  was  cloudy,  we  allowed 
satellite  matchups  from  a  day  either  side  of  our  actual 


measurements.  Another  issue  in  coastal  regions  is  that  high 
spatial  variability  can  occur  at  smaller  spatial  scales  than  are 
resolved  by  the  satellite  [Yuan  et  al. ,  2005],  further  contrib¬ 
uting  to  errors  in  satellite-based  algorithms.  Despite  poten¬ 
tial  problems,  a  great  deal  of  effort  is  now  focused  on 
advancing  satellite  applications  to  continental  shelves,  given 
that  these  are  regions  with  the  world’s  highest  oceanic 
productivity  and  which  are  increasingly  impacted  by  human 
activities.  To  obtain  the  best  satellite  retrievals  possible,  we 
have  used  a  large  in  situ  data  set  for  satellite  validation  and 
have  applied  the  most  up-to-date  atmospheric  and  absorbing 
aerosol  corrections. 

[27]  We  developed  a  multivariate  model  for  phytoplank¬ 
ton  absorption  on  the  Louisiana  shelf  from  a  suite  of 
chemical  and  physical  forcing  mechanisms,  based  on 
3  years  of  data.  The  final  multiple  regression  model 
included  13  environmental  variables  (Table  1)  and  per¬ 
formed  well  in  determining  temporal  and  spatial  variability 
in  aph  for  2002-2004  (average  r  =  0.76).  The  six  most 
important  variables,  as  determined  from  stepwise  analysis, 
were  river  discharge  with  no  lag  and  one  month  lag,  riverine 
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Figure  8.  For  April,  multiple  linear  regression  model  performance  comparing  satellite-derived  aph  (left- 
hand  column)  and  modeled  aph  (right-hand  column)  in  each  year  for  which  the  model  was  developed 
(2002-2004).  On  each  plot,  the  location  of  boxes  A  and  B  (Figure  1)  in  the  Mississippi  and  Atchafalaya 
River  plumes,  respectively,  is  indicated. 


[N03]  (no  lag),  wind  magnitude,  SST,  and  solar  radiation. 
Our  multivariate  model  is  a  statistical  model  which  allows 
for  the  determination  of  relationships  between  variables.  A 
statistical  model  does  not  necessarily  represent  the  underly¬ 
ing  causal  mechanisms  in  the  way  that  a  mechanistic  model 
would.  In  addition,  it  will  only  describe  temporal  variability 
for  the  years  of  our  study  (2002-2004)  and  on  scales 
comparable  to  the  data  scales  (monthly  in  this  case;  e.g., 
tidal  effects  would  be  difficult  to  distinguish).  For  example, 
higher  prediction  errors  would  be  observed  in  years  with 
significantly  different  forcing  variables  than  were  encom¬ 
passed  by  the  model,  as  evidenced  in  our  predictions  for 
2005,  a  year  with  river  discharge  outside  the  extremes  of 
model  parameterization  (Figure  3b).  However,  mechanistic 
models  of  phytoplankton  growth  in  river-dominated  ocean 
margins  are  still  in  relatively  early  stages  of  development,  in 
terms  of  spatial  and  temporal  resolution  [e.g.,  Green  et  al. , 
2008].  In  our  current  statistical  analysis,  we  interpret 
relationships  between  aph  and  environmental  variables 
carefully,  since  one  environmental  variable  may  be  cross- 
correlated  with  another  variable  in  our  study  (Table  2),  or 
there  may  be  important  controlling  variables  that  were  not 
considered  here. 

4.1.  Physical  Factors  Controlling  Phytoplankton 

[2s]  Freshwater  discharge  from  the  Mississippi  River  was 
one  of  the  most  important  predictors  of  phytoplankton 
absorption  over  the  hypoxic  zone.  We  emphasize  that  this 
discharge  is  related  to  buoyancy  flux  and  stratification 
effects,  and  not  with  nutrients  in  the  discharge,  because 


both  variables  are  included  separately  in  our  analysis.  For 
recent  years,  our  work  strongly  supports  “discharge  driven” 
as  one  of  the  best  descriptors  of  shelf  phytoplankton 
biomass  (Figures  5a  and  6),  as  proposed  by  Biggs  and 
Sanchez  [1997].  In  previous  work,  analysis  of  suspended 
particulate  matter  (SPM)  distributions  also  showed  that  river 
discharge  was  the  most  important  factor  determining  the 
Mississippi  River  plume’s  westward  extent  on  the  shelf 
[  Walker ,  1996].  Only  two  analyses  of  long-term  in  situ  data 
sets  have  studied  correlations  between  Mississippi  River 
discharge  and  shelf  phytoplankton  growth.  Using  data  from 
midsummer,  hypoxia  monitoring  cruises  (1985-1991),  a 
significant  positive  correlation  was  demonstrated  between 
primary  production  and  river  discharge  lagged  by  one 
month  (r  =  0.84)  at  station  C6,  a  site  located  between  the 
Mississippi  and  Atchafalaya  River  outflows  [Jus tic  et  al ., 
1993].  This  result  agrees  with  our  findings,  with  station  C6 
located  centrally  in  the  discharge-dominated  region  of  our 
study  (Figure  6).  In  a  second  study,  in  situ  primary  produc¬ 
tivity  data  from  several  years  (1988-1994)  was  averaged 
over  large  spatial  areas  of  the  shelf  [Lohrenz  et  al 1997]. 
These  average  data  sets  showed  no  significant  correlations 
with  discharge.  However,  the  contrasting  result  with  our 
study  is  likely  a  difference  in  spatial  representation,  with 
over  half  of  Lohrenz  et  al.’s  measurements  collected  outside 
of  the  discharge-dominated  regions  delineated  in  our  study. 
We  propose  two  physical  mechanisms  by  which  discharge 
could  increase  phytoplankton  biomass  on  the  shelf:  (1) 
physical  advection  of  fresh,  high-biomass  waters  farther 
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Figure  9.  For  2005,  comparison  of  (top)  satellite-derived  aph  and  (middle)  model-predicted  aph. 
(bottom)  The  %  error  between  SeaWiFS  and  modeled  images  was  calculated  using  the  absolute  value  of 
the  difference  and  referenced  to  SeaWiFS.  These  months  exemplify  application  of  the  model  to  predicting 
aph  in  a  year  that  was  not  included  in  model  development.  White  pixels  in  modeled  fields  represent 
regions  of  model  failure  (i.e.,  negative  values  of  modeled  aph ). 


down-shelf  from  river  mouths  (at  mid  to  high  salinities),  and 
(2)  increased  water  column  stratification  leading  to  reduced 
mixing  and  subsequently,  reduced  light  limitation  of  growth 
(at  low  salinities).  On  the  basis  of  previous  observations, 
light  limitation  of  phytoplankton  growth  in  the  Mississippi 
River  plume  is  probably  mainly  a  factor  at  lower  salinities 
[e.g.,  Lohrenz  et  al. ,  1999].  In  the  Louisiana  Bight,  spring/ 
summer  maxima  in  Mississippi  River  [N03]  should  rein¬ 
force  the  positive  impact  of  discharge  on  phytoplankton 
biomass  (Table  2). 

[29]  Wind  speed  was  also  an  important  contributor  to  aph 
variability  on  the  Louisiana  shelf.  The  wind-dominated 
region  was  predominantly  in  shallow  areas,  and  spatially 
distinct  from  that  of  the  discharge-dominated  region  which 
was  slightly  farther  offshore  (Figures  5  and  6).  Similarly,  a 
time  series  analysis  of  satellite-derived  SPM  concentrations 
demonstrated  spatially  distinct  zones  of  discharge-  and 
wind-driven  forcing  on  the  shelf  [Salisbury  et  al .,  2004]. 
Further  comparison  with  this  study  exemplifies  differences 
between  the  physical  forcing  of  phytoplankton  cells  that  are 
actively  growing  and  generally  less  dense,  in  comparison  to 
mineral  and  detrital  particles.  For  example,  in  our  study  of 
phytoplankton,  the  discharge-dominated  region  extended 
farther  west  from  the  Mississippi  River  mouth,  most  likely 
because  riverbome  SPM  sinks  quickly  from  surface  waters 
[e.g.,  Corbett  et  al .,  2004].  In  addition,  the  size  of  the  wind- 
dominated  region  was  much  larger  in  the  SPM  study. 
Salisbury  et  al.  [2004]  hypothesized  that  positive  wind- 


SPM  correlations  were  largely  an  expression  of  wind-driven 
processes  of  sediment  resuspension  and  transport  of  sus¬ 
pended  particles  from  shallower  to  deeper  waters.  In  con¬ 
trast,  we  observed  a  negative  correlation  between  wind 
speed  and  a We  hypothesize  that  increased  wind  speed 
causes  light  limitation  of  phytoplankton  growth  either  due 
to  a  deeper  mixed  layer  or  resuspension  of  light-attenuating 
materials  into  the  water  column.  In  fact,  changes  in  mea¬ 
sured  aph  in  shallow  regions  of  the  shelf  have  been  linked  to 
wind  events.  Following  passage  of  a  cold  front  with 
increased  wind  speeds,  aph  in  a  nearshore  region  of  the 
Louisiana  Bight  decreased  by  twofold  in  March  2002,  while 
absorption  by  CDOM  increased  [D  Sa  et  al .,  2006].  This 
change  in  aph  may  have  been  caused  by  either  the  change  in 
wind  speed  or  wind  direction  with  frontal  passage.  Because 
higher  wind  speeds  were  correlated  with  more  southerly 
winds  in  our  study  (Table  1),  the  two  factors  may  act 
together  to  decrease  aph  in  shallow  waters. 

[30]  The  importance  of  sea  surface  temperature  and  solar 
radiation  in  the  southern  region  of  the  hypoxic  zone 
indicates  a  contribution  from  upwelled  nutrients  to  phyto¬ 
plankton  farther  offshore  on  the  Louisiana  shelf.  Whereas 
nutrient  concentrations  are  higher  in  low-salinity  shelf 
waters,  farther  offshore  on  the  shelf  nutrient  concentrations 
are  often  negligible,  especially  in  Fall  and  Winter  [e.g., 
Lohrenz  et  al .,  1999;  Chen  et  al .,  2000].  In  our  study,  both 
solar  radiation  and  SST  varied  seasonally  with  minima  in 
winter.  In  contrast,  in  the  southern  region  of  the  hypoxic 
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Month  in  2005 

Figure  10.  For  2005,  time  series  comparisons  of  satellite- 
derived  and  model-predicted  aph  in  the  (a)  Mississippi  River 
plume  (Box  A)  and  (b)  Atchafalaya  River  plume  (Box  B). 
Both  satellite  and  modeled  time  series  peaked  in  spring  to 
early  summer  in  the  Mississippi  plume  and  in  midsummer 
in  the  Atchafalaya  plume. 


zone,  aph  often  displayed  wintertime  peaks,  and  was  thus 
negatively  correlated  with  the  two  forcing  variables.  In  deep 
Gulf  of  Mexico  waters,  previous  studies  have  shown  that 
thermally  driven  mixed  layer  depths  were  the  main  factor 
driving  the  winter  maxima  in  phytoplankton  biomass 
[ Walsh  et  al.,  1989;  Miiller-Karger  et  al.,  1991].  The 
cooccurrence  of  high  algal  biomass,  cooler  SST’s,  and 
deeper  mixed  layers  provided  strong  evidence  that  primary 
production  in  the  open  Gulf  was  controlled  by  the  upward 
flux  of  N03.  We  hypothesize  that  the  same  thermal  mech¬ 
anism  which  controls  phytoplankton  growth  in  the  open 
Gulf  also  plays  a  role  in  offshore  shelf  waters.  Previously, 
studies  have  shown  upwelling  favorable  currents  to  induce 
bottom  Ekman  upwelling  that  injects  nutrient-rich  deep 
water  onto  the  shelf,  causing  localized  high  surface  values 
of  nutrients  and  phytoplankton  biomass  [Fox  et  al.y  1987; 
Lopez-Veneroni  and  Cifuentes ,  1994;  Sahl  et  al .,  1993; 
Chen  et  al.y  2000].  Similarly,  for  our  study  years  (2002- 
2004),  we  observed  wintertime  winds  and  currents  that 
caused  circulation  to  flow  off-shelf,  providing  upwelling 
favorable  conditions.  On  the  basis  of  our  study,  both  deeper 
mixed  layers  and  favorable  circulation  patterns  in  winter 
appear  to  promote  the  vertical  entrainment  of  nutrients  from 
either  benthic  or  upwelled  sources  and  to  increase  phyto¬ 
plankton  biomass  over  the  hypoxic  zone.  Thus,  we  identify 
two  different  mixing  regimes  on  the  shelf:  (1)  in  nearer 
shore  waters,  mixing  negatively  affected  phytoplankton 
biomass  due  to  increased  light  limitation  (Figure  6,  dis¬ 
charge-  and  wind-dominated  regions),  and  (2)  farther  off¬ 
shore,  mixing  positively  affected  biomass  due  to  vertical 
entrainment  of  nutrients  (Figure  6,  radiation-  and  SST- 
dominated  regions).  Similarly,  multiple  sources  of  nutrients 
to  phytoplankton  (i.e.,  riverine  and  upwelled)  have  been 
documented  in  other  river-dominated  margins,  such  as  the 
Amazon  River  plume  [e.g.,  DeMaster  and  Pope ,  1996]. 


4.2.  Chemical  Factors  Controlling  Phytoplankton 

[31]  Over  the  hypoxic  zone,  Mississippi  River  nitrate 
concentration  was  one  of  the  most  important  predictors  of 
phytoplankton  absorption.  Riverine  [NO3]  was  most  impor¬ 
tant  in  nearer  shore  waters  (Figures  5c  and  6).  As  in  our 
above  analysis  of  discharge,  only  two  analyses  of  long-term 
in  situ  data  sets  have  studied  correlations  between  Missis¬ 
sippi  River  NO3  and  phytoplankton  growth  in  Louisiana 
coastal  waters.  For  an  eight  year  period  (1985-1992), 
measurements  at  Station  C6  on  the  Louisiana  shelf  (station 
location  shown  on  Figure  6)  showed  that  primary  produc¬ 
tivity  in  the  summertime  was  significantly  correlated  with 
riverine  N03  flux  lagged  by  one  month  (r  =  0.73  [Justic 
et  al.y  1997]).  However,  the  authors’  previous  research  for 
approximately  the  same  years  (1985-1991)  and  same 
location  resulted  in  a  higher  correlation  between  primary 
productivity  and  discharge  lagged  by  one  month  (r  =  0.84 
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Figure  11.  Spring  model  experiment  showing  variability 
in  aph  caused  by  changing  all  environmental  variables, 
except  N03  flux,  with  comparison  of  springtime  (a)  average 
aphy  (b)  maximum  aph ,  and  (c)  minimum  aph.  Average  aph 
was  determined  by  inputting  average  springtime  (April— 
June,  2002-2004)  environmental  variables  into  the  model. 
Maximum  and  minimum  aph  were  modeled  by  using  the 
limits  of  all  environmental  variables,  except  NO3  flux, 
during  this  time  period.  The  maximum  and  minimum 
scenarios  were  determined  on  the  basis  of  an  average  aph 
calculated  over  the  entire  region.  Thus,  aph  may  not  be 
higher  in  the  maximum  scenario  at  specific  locations,  such 
as  between  Atchafalaya  and  Terrebonne  Bays. 
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[Justic  et  al. ,  1993]).  Hence,  we  conclude  that  river  dis¬ 
charge  was  more  important  than  N03  concentration  in 
driving  productivity  at  station  C6,  in  concurrence  with  our 
findings,  as  presented  above.  In  contrast,  Lohrenz  et  aV s 
[1997]  study  around  the  delta  showed  that  both  riverine 
N03  flux  and  N03  concentration  were  important  predictors 
of  primary  productivity  during  a  7-year  period  ( 1 988  — 
1994).  Many  of  their  measurements  were  in  a  region  just 
southwest  of  the  delta  where  [N03]  was  often  also  an 
important  predictor  of  aph  in  our  study  (Figure  6).  The 
generally  good  agreement  between  the  results  of  our  study 
and  the  limited  correlations  presented  in  previous  field  work 
is  somewhat  surprising.  It  can  be  difficult  to  directly 
compare  results  derived  from  point  measurements  to  those 
derived  from  more  synoptic  satellite  imagery,  especially 
given  differences  in  time  period  between  our  study  and 
these  previous  measurements.  As  well,  we  are  examining 
variability  in  aph>  which  is  more  an  indicator  of  phytoplank¬ 
ton  biomass  than  a  rate  process,  such  as  primary  production. 
Though,  phytoplankton  biomass  and  primary  production 
tend  to  be  strongly  correlated  in  this  region  [Lohrenz  et  aly 
1999;  Chen  et  al ,  2000]. 

[32]  One  of  our  main  findings  was  that  riverine  nitrate 
concentration  and  river  discharge  were  equally  important 
determinants  of  aph  on  the  shelf  (Figure  6).  Previous 
observations  have  demonstrated  that  riverine  effects  on 
salinity  (and  thus  buoyancy  and  stratification)  are  wide¬ 
spread  along  the  shelf  [Dinnel  and  Wiseman ,  1986],  where¬ 
as  riverine  effects  on  nutrient  distributions  are  confined 
closer  to  the  river  mouths  [e.g.,  Sahl  et  al.y  1993].  In  fact, 
we  showed  here  that  both  of  these  chemical  and  physical 
processes  are  very  important  in  determining  surface  phyto¬ 
plankton  distributions,  and  that  their  effects  are  observed  in 
unique  regions  of  the  shelf  -  generally  in  nearer  shore 
waters  for  [N03],  compared  to  river  discharge.  Possible 
relationships  with  the  in  situ  production  of  ammonium  may 
help  explain  the  importance  of  riverine  [N03]  away  from 
river  mouths,  such  as  in  the  Atchafalaya  River  plume 
(Figures  5c  and  6).  On  the  basis  of  the  correlation  and 
regression  analyses  presented  here,  we  also  suggest  that  the 
vertical  entrainment  of  nutrients  from  benthic  regeneration 
and/or  upwelled  sources  into  surface  shelf  waters  is  region¬ 
ally  important  to  phytoplankton  growth  (e.g.,  in  the  more 
offshore  region  of  the  hypoxic  zone).  In  the  case  of  benthic 
inputs,  the  Mississippi  and  Atchafalaya  Rivers  may  have 
been  the  initial  source  of  these  nutrients. 

[33]  We  observed  limited  spatial  influence  of  silicate  and 
phosphate  concentrations  on  phytoplankton  absorption  on 
the  Louisiana  shelf.  The  role  of  [P04]  is  somewhat  less  clear 
than  that  of  [Si03],  because  [P04]  was  highly  correlated 
with  SST,  whereas  [Si03]  generally  had  low  correlations 
with  other  environmental  forcing  variables  (Table  2).  Our 
final  multivariate  model  did  not  include  [PO4],  but  did 
include  [Si03]  for  its  role  in  forcing  aph  west  of  Atchafalaya 
Bay.  Although  nitrogen  often  controls  primary  production 
in  shelf  waters  near  the  Mississippi  River  delta,  phosphorus 
and  silicate  limitation  have  also  been  observed  [e.g.,  Dortch 
and  Whitledge ,  1992;  Lohrenz  et  al.y  1999].  In  this  study,  we 
have  only  considered  nutrient  inputs  from  the  Mississippi 
River,  whereas  nutrient  contributions  from  other  sources 
may  lead  to  a  somewhat  different  understanding  of  nutrient 
effects  on  phytoplankton.  For  example,  processes  of  bio¬ 


logical  uptake,  remineralization  of  organic  matter,  and/or 
upwelling  may  all  affect  in  situ  nutrient  concentrations  and 
which  nutrients  are  limiting  to  growth  [e.g.,  Lohrenz  et  al.y 
1999].  However,  our  results  should  be  fairly  robust  in 
nearshore  waters  that  are  directly  impacted  by  riverine 
outflow.  In  these  waters,  phytoplankton  biomass  was  the 
highest  on  the  shelf  and  was  clearly  correlated  with  riverine 
[N03]  (Figure  5c),  but  to  only  a  much  lesser  extent  with 
either  riverine  [P04]  or  [Si03].  Thus,  we  interpret  our 
results  to  indicate  that  riverine  [N03]  is  the  primary  nutrient 
controlling  phytoplankton  variability  on  the  Louisiana  shelf, 
in  concurrence  with  previous  findings  [Turner  et  a/.,  2006]. 

43.  Implications  for  Management  of  Hypoxia 

[34]  Stratification  and  water  column  mixing  play  an 
important  role  in  determining  phytoplankton  biomass  and 
the  areal  extent  of  hypoxia  on  the  Louisiana  shelf.  Seasonal 
hypoxia  is  the  combined  result  of  vertical  stratification  of 
the  water  column  and  respiration  of  surface  organic  carbon 
in  bottom  waters.  Much  attention  has  focused  on  eutrophi¬ 
cation  by  increased  N03  loading  over  the  last  50  years  as 
the  main  cause  of  hypoxia  in  Louisiana-Texas  shelf  bottom 
waters  [e.g.,  Turner  et  al.y  2006,  and  references  therein]. 
However,  it  is  clear  that  physical  mechanisms  also  strongly 
control  the  extent  of  hypoxia  and  that  [N03]  is  not  the  sole 
cause.  Previous  work  described  a  large  region  of  the  shelf 
which  is  affected  by  freshwater  flow,  where  the  euphotic 
zone  is  nitrogen  limited  and  stratification  plays  the  most 
important  role  in  determining  hypoxia  [Rowe  and 
Chapman ,  2002].  The  authors  emphasized  that  stratifica¬ 
tion-induced  hypoxia  has  been  observed  in  other  coastal 
areas  (e.g.,  the  Pamlico  River  Estuary  [Stanley  and  Nixon , 
1992])  and  that  Louisiana  shelf  hypoxia  correlates  just  as 
well  with  the  strength  of  the  density  gradient  as  it  does  with 
[N03].  Recent  results  from  the  combined  use  of  a  hydro- 
dynamic  and  oxygen  respiration  model  concluded  that 
hypoxia  formation  is  primarily  dependent  on  vertical  mix¬ 
ing  processes,  as  determined  by  freshwater  input,  bottom 
topography,  and  winds  [Hetland  and  DiMarco ,  2008]. 

[35]  Our  findings  further  support  the  importance  of 
mixing  on  the  shelf.  Although  we  observed  spatial  hetero¬ 
geneity  in  the  physical  variables  driving  phytoplankton 
biomass  on  the  shelf  (Figure  6),  a  common  factor  among 
all  of  these  physical  variables  was  that  their  impact  on 
phytoplankton  is  likely  a  result  of  changes  in  stratification. 
For  example,  correlation  analyses  showed  that  in  nearshore 
and  midshelf  waters,  lower  wind  speeds  and  increased 
freshwater  discharge  both  increased  phytoplankton  bio¬ 
mass  (e.g.,  Table  1).  In  offshore  shelf  waters,  a  thermally 
driven  mixing  regime  was  observed  where  we  hypothesize 
that  increased  mixing  in  winter  positively  affected  phyto¬ 
plankton  biomass  due  to  the  vertical  entrainment  of 
nutrients  (Figure  6).  As  well,  the  obvious  impact  of  mixing 
is  observed  in  the  passage  of  storms  over  the  shelf.  For 
example,  in  the  years  of  our  study,  the  areal  extent  of 
hypoxia  was  significantly  reduced  in  years  when  tropical 
storms  and  hurricanes  had  a  stronger  impact  on  our  study 
area  (2003  and  2005,  <12,000  km2)  compared  to  years 
with  lower  storm  activity  (2002  and  2004,  >15,000  km2). 
We  also  observed  lower  aph  in  July  of  2003  and  2005,  the 
years  with  higher  wind  speeds  (Figures  3a  and  3c). 
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[36]  A  significant  degree  of  variability  in  phytoplankton 
biomass  on  the  Louisiana  shelf  can  be  expected  even  for 
constant  riverine  N03  flux.  A  springtime  experiment  with 
our  model  suggested  that  average  aph  over  the  hypoxic  zone 
could  vary  by  an  average  37%  owing  to  changes  in 
environmental  variables  other  than  N03  flux,  and  by  more 
in  some  regions  of  the  shelf  (e.g.,  LA  Bight;  Figure  11). 
Understanding  the  various  factors  which  influence  hypoxia 
is  especially  important  given  similarity  in  the  natural 
variability  we  observed  in  aph  (37%)  and  the  proposed 
30%  reduction  in  annual  nitrogen  flux  into  the  Gulf  of 
Mexico,  as  recommended  by  the  Federal-State-Tribal  Ac¬ 
tion  Plan  for  reducing,  mitigating,  and  controlling  hypoxia 
[Mississippi  River/Gulf  of  Mexico  Watershed  Nutrient  Task 
Force ,  2001].  The  management  and  economic  implications 
for  reaching  this  nitrogen  reduction  goal  are  far-reaching, 
involving  potential  large-scale  changes  in  agricultural  fer¬ 
tilization  practices  in  Midwestern  states  [Mitsch  et  al., 
2001].  Hence,  it  is  important  that  scientists  and  environ¬ 
mental  managers  can  explain  the  resulting  changes  in 
hypoxia  size  as  nutrient  flux  is  decreased.  To  the  extent 
that  there  is  a  linkage  between  surface  production  and 
hypoxia,  our  spring  model  experiment  suggests  that  changes 
in  hypoxia  size  can  be  expected  from  variability  in  certain 
physical  factors.  For  example,  physical  factors  likely  con¬ 
tributed  to  the  2002  hypoxic  zone  being  the  largest  in 
history.  In  spring  2002,  a  year  with  only  moderate  spring- 
summer  N03  concentrations  (Figure  3d),  we  observed 
stronger  winds  and  currents  directed  to  the  west,  compared 
to  either  2003  or  2004,  which  caused  retention  of  freshwater 
and  phytoplankton  biomass  on  the  shelf  (Figure  8).  In 
predictive  models  of  hypoxia  size,  there  is  currently  little 
consensus  as  to  which  physical  processes  are  most  impor¬ 
tant  [e.g.,  Scavia  et  al .,  2004].  We  suggest  that  the  predic¬ 
tive  capability  of  these  models  would  be  improved  by 
incorporating  some  of  the  physical  factors  that  we  have 
identified  as  important  on  the  shelf,  especially  that  of  water 
column  mixing. 

5.  Conclusions 

[37]  Variability  in  phytoplankton  biomass  in  Louisiana 
shelf  surface  waters  is  driven  by  a  host  of  chemical  and 
physical  forcing  variables.  Of  these  factors,  we  conclude 
that  recent  variability  in  satellite-derived  aph  (2002-2004) 
is  most  strongly  driven  by  Mississippi  River  discharge  and 
riverine  [N03],  followed  by  SST,  solar  radiation,  and  wind 
speed.  On  the  basis  of  these  findings,  two  different  mixing 
regimes  were  identified  on  the  shelf,  one  in  nearer  shore 
waters  where  mixing  negatively  affected  phytoplankton 
biomass  due  to  light  limitation  and  a  second  region  farther 
offshore  where  mixing  positively  affected  biomass  due  to 
upwelled  nutrients.  This  offshore  region  indicates  vertical 
entrainment  as  a  source  of  nutrients  to  shelf  phytoplankton, 
though  likely  of  secondary  importance  to  riverine  nutrient 
inputs  [e.g.,  Biggs  and  Sanchez ,  1997].  A  multiple  regres¬ 
sion  model  was  developed  from  13  of  the  environmental 
variables  initially  included  in  our  analysis,  and  performed 
well  in  modeling  spatial  and  temporal  aph  distributions  in 
2002-2004,  as  well  as  in  predicting  aph  in  2005.  We 
propose  that  further  insights  into  controls  on  shelf  phyto¬ 
plankton  biomass  could  be  gained  from  development  of 


seasonal  models,  as  have  been  used  in  predicting  suspended 
sediment  distributions  on  the  shelf  [e.g..  Walker ,  1996].  As 
well,  given  inherent  difficulties  with  satellite  algorithms  in 
coastal  regions,  satellite  retrievals  of  phytoplankton  biomass 
in  this  region  may  be  further  improved  through  the  use  of 
semianalytical  algorithms  [e.g.,  Lee  et  al .,  2002;  Maritorena 
et  al .,  2002],  and  continued  validation  with  in  situ  obser¬ 
vations.  Finally,  we  found  that  the  inclusion  of  currents 
from  a  Gulf  of  Mexico  circulation  model  improved  our  final 
model  for  aph.  We  suggest  that  future  research  will  benefit 
from  improved  circulation  models  for  river-dominated  shelf 
regions,  including  higher  spatial  resolution,  incorporation  of 
actual  (rather  than  climatological)  river  inputs,  and  better 
parameterization  of  surface  salinity  flux.  In  summary,  our 
findings  show  that  time  series  analyses  of  ocean  color 
imagery  can  provide  valuable  insights  into  phytoplankton 
dynamics  in  shallow  shelf  waters  of  the  Gulf  of  Mexico.  As 
well,  we  have  demonstrated  that  remote  sensing  in  combi¬ 
nation  with  environmental  forcing  mechanisms  may  provide 
scientists  and  environmental  managers  with  powerful  tools 
for  understanding  future  changes  in  hypoxia  on  the  Louisiana 
shelf. 
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